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PACS 04 . 20 . -q - Classical general relativity 

PACS 04 . 25 . -g - Approximation methods; equations of motion 

PACS 97 . 80 . -d - Binary and multiple stars 

Abstract. - Based on the Arnowitt-Deser-Misner (ADM) canonical formulation of general rel- 
ativity, a canonical formulation of gravitationally interacting classical spinning-object systems is 
given to linear order in spin. The constructed position, linear momentum and spin variables ful- 
fill standard Poisson bracket relations. A spatially symmetric time gauge for the tetrad field is 
introduced. The achieved formulation is of fully reduced form without unresolved constraints, 
supplementary, gauge, or coordinate conditions. The canonical field momentum is not related to 
the extrinsic curvature of spacelike hypersurfaces in standard ADM form. A new reduction of the 
tetrad degrees of freedom to the Einstein form of the metric field is suggested. 
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The canonical formulation of spinning objects under 
gravitational interaction is an important issue in general 
relativity. For the gravitationally interacting Dirac field, 
e.g., several investigations have been undertaken to set- 

: tie the problem, 1-6 . For spinning classical objects like 
rotating black holes, neutron stars, or other stars simi- 

: lar advances have not yet been obtained. For objects in 
external gravitational fields, canonical formulations were 
constructed by Kiinzle [7J and, very recently, by Barausse 
et al. [8]. In the paper by Yee and Bander [5] a Routhian 
approach was developed with the spin part in canonical 
form. Within effective field theory techniques, the lat- 
ter approach was pushed forward by Porto and Rothstein 
[10] to self-gravitating spinning objects operating in non- 
reduced phase space of the matter variables with the re- 
duction to be performed within post-Newtonian (PN) ap- 
proximations. The formalism was shown to operate well to 
the next-to-leading order spin-spin coupling. Based on the 
Arnowitt-Deser-Misner (ADM) canonical formulation of 
general relativity, , starting from the stress-energy ten- 
sor for pole-dipole objects by Tulczyjew and Dixon, a fully 
reduced matter-field canonical formalism was achieved, 
but to next-to-leading order in the PN approximation 
only [12]. Most recently, this approach has been gener- 
alized to a next-order PN approximation [13] . Within the 
previously developed canonical scheme even quadratic-in- 
spin interactions could be treated successfully, [14] , The 
leading order spin-gravity interaction has found various 
derivations, for quantum and classical spinning particles 
see, e.g., [15] . for black holes and other bodies, [16]. With 



the aid of a test-spin-type approach, the next-to-leading 
order spin-orbit dynamics in canonical form has been ob- 
tained for the first time only recently, [17]. It has turned 
out fully equivalent to the non-canonical result in [T8] . 

In the present article, a canonical formulation for self- 
gravitating spinning classical objects will be developed 
which, to linear order in the spin, is valid for arbitrary 
gravitational fields. The formulation makes crucial use of 
a spatially symmetric time gauge for the tetrad field sim- 
ilar to the one introduced by Kibble [2] and, like in Kib- 
ble's investigation of the gravitationally interacting Dirac 
field, a full reduction is achieved without remaining un- 
resolved constraints. The gravitational field turns out to 
be treated within a field-momentum-generalized version 
of the ADM canonical formalism. Thus, about 50 years 
after the seminal work by ADM, classical spinning-object 
systems receive a canonical implementation into the Ein- 
stein theory of gravity. The importance of a canonical 
formulation is discussed, e.g., in [19]. Most notably it 
allows a thorough analysis of the consistency of the sys- 
tem of equations consisting of the Mathisson-Papapetrou 
equations, ^ and (J7J), and the Einstein equations with 
the Tulczyjew-Dixon stress-energy tensor (0) as a source. 
(This system is applicable to the dynamics of compact ob- 
jects like black holes or neutron stars if rotation is not too 
rapid, i.e., if spin-squared terms can be neglected, and if 
tidal deformation has no effect. It was thus the basis for 
many investigations of such objects in the past. Of course, 
also weakly gravitating objects can be modeled in this way 
as long as their deformation is negligible.) In particular, 
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the action used in the present paper belongs to the class 
of so-called derivative coupled theories, which are known 
to potentially reveal subtle inconsistencies when formu- 
lated canonically, see [TH]. However, for the Dirac field 
the problematic terms containing the extrinsic curvature 
of the space-like hypersurfaces can be eliminated by a re- 
definition of the Dirac field. In our case, those problematic 
terms are absorbed into the canonical momentum of the 
matter. Besides rather mathematical consistency consid- 
erations that now become available, applications of the 
canonical formulation given in the present paper to the 
post-Minkowskian (see, e.g., [50]) and to higher order PN 
approximations (next-to-next-to-leading order spin-obit in 
particular) are most interesting. Further, the suggested 
gravitational-field reduction will apply to the Dirac field 
as well with possibly nice prospects for applications, cf., 

mi- 

Latin indices from the middle of the alphabet are run- 
ning through i = 1,2,3. We utilize three different frames 
here, denoted by different indices. Greek indices refer to 
the coordinate frame and have the values /i = 0,i. Lower 
case Latin indices from the beginning of the alphabet re- 
fer to the local Lorentz frame, while upper case ones de- 
note the so called body-fixed Lorentz frame. The values 
of these Lorentz indices are marked by round and square 
brackets as a = (0), (i) and A — [0], [i], respectively, e.g., 
A = [0], [1], [2], [3]. Partial derivatives are denoted ilt . 
Both the speed of light, c, and the Newton gravitational 
constant, G, are put equal to 1. 

Our starting point will be an action functional W 
which is invariant against general four-dimensional coor- 
dinate transformations, general local tetrad rotations, and 
reparametrization of the objects affine time parameter. In 
terms of a Lagrangian density C the action reads 

W[e alx ,z»,p ll ,h Ca ,S ab ,\ a 1 ,\ m M = Jd 4 xC, (1) 

and must be varied with respect to the tetrad field e a/i , 
the Lagrange multipliers A", X^w, A3, position z M and lin- 
ear momentum p^ of the object, as well as with respect 
to angle-type variables A Ca and spin tensor S ab associ- 
ated with the object. The angle variables are represented 
by a Lorentz matrix satisfying A Aa A Bb r]AB — V ab or 
A Aa A Bb ?f b = r] AB , where r] AB = diag(-l, 1, 1, 1) = i] ab , 
which must be respected upon variation, see [22]. The 
matter part of the Lagrangian density reads 



C 



M 



dr 



1 c , , ab 



dz^ 
"dr" 



1 d9 ab 



(4) , 

(2) 

with dd ab = A c a dA cb = -d9 ba and <5 (4) = 8{x v - z v {t)) 
the 4-dim. delta function, J d 4 x<5(4) = 1. The objects 
affine time variable is r. The Ricci rotation coefficients 



are given by uJ m p = e aa e b pbj^ 



= -T 



with rfl 

fel symbols of first kind and 



(4) 



Kfl'/Sa.M +9PiJ,,a " 9an,p) the 4-dim. Christof- 



ebuV 



ab 



the 4-dim. 



metric. As in [22], the matrix A Ca can be subjected to 
right (or left) Lorentz transformations, which correspond 
to transformations of the reference frame (or the body- 
fixed frame) . The spin part of Cm can be obtained from 
the kinetic term \S ab ^r in Minkowski space, see [22], by 
promoting the global symmetry of this term under right 
Lorentz transformations (i.e., transformations of the ref- 
erence frame), to a local symmetry, a usual procedure for 
the construction of gauge theories. This is achieved by 
introducing the ui^ lb term in ([2]). The matter constraints 
are given by 



c 



dr 



X a lP b S ab + \ m A® a p a ~^-( P 2 + m 2 ) 



?(4) 



where m is the constant mass of the object and p 2 = p^p^ ■ 
The constraint p b S ab = (spin supplementary condi- 
tion, SSC) states that in the rest frame the spin ten- 
sor contains the 3-dim. spin Su\(j<\ only (i.e., the mass 
dipole part S(o)(i) vanishes), while the conjugate con- 
straint AM°p = ensures that A Ca is a pure 3-dim. ro- 
tation matrix in the rest frame (no Lorentz boosts), see 
[22] , Finally, the gravitational part is given by the usual 
second order Einstein-Hilbert Lagrangian density 



(4) 



(TD), 



(4) 



where g is the determinant of the 4-dim. metric and Pi/ 4 ^ 
is the 4-dim. Ricci scalar. Using a second order form of 
the gravitational action, i.e., not varying the connection 
independently, ensures that the torsion tensor must van- 
ish, see, e.g., [5|. A total divergence in the Lagrangian 
density may be added without affecting the equations of 
motions, as they can be obtained from local variations. 
The term denoted (TD) is written here as a reminder of 
this fact, and will become important later on. The com- 
plete Lagrangian density is the sum 
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C 



M 



(5) 



We assume asymptotical flatness as a boundary condition 
of the spacetime. 

Variation of the action SW = leads to the equations 
of motion (EOM) for the matter variables 



D5 ab Q DA Ca 



0, it* 



dz^ 



Dt Dr dr 

Dr ~ 2 ' 



A 3 P M - 



(6) 
(7) 



as well as to the usual Einstein equations with the stress- 
energy tensor 



e£ 5(Cm + C c ) 



Se a 



= / dr 



9/\\a 



(8) 
, (9) 
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where R^) ab is the 4-dim. Riemann tensor, 



| a 



denotes the 



4-dim. covariant derivative, and d and D, respectively, 
denote ordinary and covariant total derivatives. Here it 
was already used that preservation of the constraints in 
time requires A" to be proportional to p a and A 2 [i] to be 
zero, so that A" and X2U] drop out of the matter EOM 
and the stress-energy tensor. The Lagrangian multiplier 
A3 = Aa(r) represents the reparametrization invariance of 
the action (notice A3 = y— u 2 /m). Further, an antisym- 
metric part of the stress-enery tensor 



2 V=5 



\\p 



= / drlB^ 
2 Dt 



fJ 



= (10) 



vanishes and it holds 



by virtue of the mat- 
ter EOM. Obviously, the spin length s defined by 2s 2 = 
S a bS ab is conserved. 

To our knowledge, the Lagrangian density in the form 
given here is used for the first time, though the expres- 
sions applied in, e.g., [3 El UHl are somewhat related. 
In [TJE] spinning objects in an external gravitational field 
are treated, in contrast to self-gravitating spinning objects 
here. In [23] the metric is the fundamental variable (e.g., 
used in functional integrations) instead of the tetrad field, 
and [TO] gives a Routhian version of [23] . Varying with 
respect to A AfJl {e J in [33]), however, is quite subtle as one 
must respect A^Aal* — g^v That is, the variations SA^ 
and Sg^v are not independent. The tetrad field e aM thus 
implicitly also appears in [23] . but it is fixed as a func- 
tional of the metric, see eq. (34) in [33J. In this paper, 
however, it is crucial that Se afl and 5A Ca are independent 
and that the full-fledged gauge freedom of the tetrad is 
manifestly available. 

The approach in this paper to a fully reduced canon- 
ical framework is to eliminate all constraints and gauge 
degrees of freedom from the action. After that, only the 
independent variables, which parametrize the constraint 
surface, must be varied. The action is then transformed 
into canonical form by certain variable transformations. 
A (3+l)-split with respect to a timelike unit 4-vector n^ 
with components = (-iV, 0,0,0) or n^ = (1, -N*)/N, 
where N is the lapse function and N l the shift vector, most 
naturally fits to a fully reduced canonical formulation of 
gravity. The three matter constraints can be solved as 



np = n^p^ = -Jm z + J^PiPj > 



(11) 



nS t = n^S^ = Pkl 3Sji = gijnS^ , (12) 



Dp 



AbK°) = aM»M a [0]q = -— 



P 



(0) 



(13) 



in terms of pi, Su\rj\, and A^ k \ so one can put Cc = 
now. Here r ) lJ is the inverse of the induced 3-dim. metric 
gij = jij of the hypersurfaces orthogonal to n M . A split of 
the Ricci rotation coefficients results in 



(14) 



n"w M = K H - 9ij ^ + ^(el k - et k N l ) , (15) 



Ni 



(16) 



n^i = KijN' - N ;i - 9i f-f + ^(eg, - e^N 1 ) , 

(17) 

where ; j denotes the 3-dim. covariant derivative, Tjn~ the 
3-dim. Christoffel symbols and the extrinsic curvature Kij 
is given by 2NKij = —gij.o + 2Nu. i j\. For convenience, we 
will immediately go to the time gauge e't-. = n* 1 , see [2"4"] . 
also [TJ[2l[5], as lapse and shift then turn into Lagrange 
multipliers, like in the ADM formalism (e.g., the N 3 k /N 

and N J /N terms are canceled in (fTB"]) and (fTT]) ). It holds 



(19) 
(20) 



In passing we mention that in [7] a completely different 
tetrad field has been chosen for a canonical formulation of 
test spinning particles moving in gravitational fields. 

The matter action in the covariant SSC p b S a b = turns 
into 

C M = Cmk + C MC + Cgk + (td) , (21) 

where (td) denotes an irrelevant total divergence. The 
terms attributed to the kinetic matter part are given by 



(0) 

e i 


= = 4), 


e (0) 


= N= l/e° 0) 


N' 


= -Ne\ 0) , 


(») 


= N*ef , 


9ij 


(m) 

— e i e (m)j , 




= elm/ m) > ■ 



Cmk = 



p. + K^nS 3 +A kl e (j)k e\ 



1 

2 np 



z 6 + - — pid 
Znp 



S, 



(i)U) 



a a -1 a(») A MO') 

nS {i) p U) - nS (j) p (i) A [fc] A 



np 



with A ZJ defined by 

gikgjiA kl 



2np 



-5, 
(22) 

(23) 



The delta function S is defined as J d 3 xS(x l — z l (t)) = 1, 
and the gauge r = z° = t was chosen, where t is the time 
coordinate of the object. A "dot" denotes the derivative 
with respect to t. The complicated structure of these ki- 
netic terms represents the Dirac bracket arising from the 
covariant SSC. The matter parts of the gravitational con- 
straints result from 



zi/mattcr 



C MC = —NTL™ 1 
with 

^matter = § _ Rl jPf^l S _ f n S k 5). k . 

np 



(24) 



(25) 
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fa + K ijn S j )S 



(1 rnkn r . r(* l)mPk r ''Sl 

-7 SikS + Sl 7 ; 5 

\2 np 



(26) 



These coincide with the densitized projections % matter = 
V^fT^n^n" andKf atter = -^T iv n v of the stress-energy 
tensor in covariant SSC, eq. ([9j, see also [12], where 7 = 
det(Tjj). Further, some terms attributed to the kinetic 
part of the gravitational field appear as 

£ GK = A»e {k)i e$6. (27) 

Now we proceed to Newton- Wigner (NW) variables z l , 
Pi, S(i)(j\, and Ay^', which turn the kinetic matter part 
Cmk into canonical form. The variable transformations 
read 



nS l 



m — np 
PinSj 



nSi 



m — np m — np 
,(3) 



(28) 
(29) 



A WW) = AWW ( 5 kj + P{k)P ' - s ) , (30) 
\ m(m — np) J 

(31) 



■-Tip) J ' 

Pi = Pl + K ijn S ] + A kl e u)k e^) 
1 

.2 k ' np 



where A lJ is given by 



A kl 

gikgjiA = -bi : 



np(m — np) 



(32) 



The NW variables have the important properties 



2s 2 = const and Af^A^Ki) 
[k] 



Sij, which 



implies that dflWW = AfSdA^^ i s antisymmetric. The 
redefinitions of position, spin tensor, and angle-type vari- 
ables are actually quite natural generalizations of the 
Minkowski space versions, cf., refs. [22] 125] . to curved 
spacetime. However, there is no difference between lin- 
ear momentum pi and canonical momentum P$ in the 
Minkowski case. In these NW variables, one has 



£ 



GK 



£ 



MK 



£. 



GK 



£ 



MK 



(td), 



with (from now on S = S(x l — z l (t))) 



jC MK = P i k i 6 + -S m) 9^5, 
£ GK =A^e {k)i e%5. 



(33) 

(34) 
(35) 



Notice that all j>i terms in the action have been canceled 
by the redefinition of the position. Further, all terms 
were eliminated from £mc an d £mk by the redefinition 
of the linear momentum. If the terms explicitly depend- 
ing on the triad ejj are neglected, the known source terms 
of Hamilton and momentum constraints in canonical vari- 
ables, respectively (4.23) and (4.25) in [T2], are obtained. 



The final step goes with the ADM action functional of 
the gravitational field, [TT][2j)][27], but in tetrad form as 
derived in [3] • In passing we mention that Kibble, [2] , was 
applying the Schwinger canonical formalism which uses a 
different set of field variables, [21] ■ The canonical momen- 
tum conjugate to &i k )j is given by 



d£ 



de 



(k)j,o 



3 j*V J ' + e\ k) 8irA ij 6, 



where 



Legendre transformation leads to 



(36) 



(37) 



-GK 



£ ( 



8n 



i" (fcW e ( fc) i>0 - 



16tt 



GC 



(td). (38) 



The explicit form of the non-irrelevant total divergence E^i 
emerges from the total divergence that can be added to 
the action, see eq. As shown in [11,26,27 by different 
methods, see also [24] . it must be given by £i — gijj — gjjj 
for asymptotically flat spacetimes. The total energy reads 
E = y-U § d 2 Si£i. The constraint part takes the form 



£ 



GC 



-NU 



field 



field 



(39) 



with 



H 



field 



1671-^/7 



7R+2 {m^f 



ik 

gijgkiv 7T 



o,field _ _L fl .._i* 

n i — 8ir ;fe ' 



(40) 
(41) 



where R is the 3-dim. Ricci scalar. Due to the symmetry of 
ir l \ not all components of 7r' fc " are independent variables 
(i.e., the Legendre map is not invertible), leading to the 
additional constraint — %ixA- % ^b. This constraint will 
be eliminated by going to the spatial symmetric gauge 
0)3 = e ij = e 3\ Then the triad 



is fixed as the matrix square-root of the 3-dim. metric 

eije jk = g lkl or 

en = \/{gki) ■ 



Therefore, we can define a quantity Bfj as 

Tiki 



or, in explicit form, 



2B^ 



de r , 



dg, 



ki 



dg, 



kl 



(42) 
(43) 

(44) 



This expression may be evaluated perturbatively, cf. 
It holds Bfj5 k i = 0. Furthermore, 



which gets applied as 



(45) 



(46) 
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with the new canonical field momentum 

<ln = *** + ^A {lj) 8 + l^B%A [kl] 5. (47) 

Notice that (|45l) can also be used to replace the triad terms 
in l|gl]l. 

Next the gravitational constraints arising from the vari- 
ations 5N and SN l , 

H md + ^matter = q ^ ^ficld + ^matter = q ^ ^ 

are eliminated by also imposing the gauge conditions 

<i,,. <»• (49) 

which allow for the decompositions 



J, 



/ / / _,{ZgPaj Z 3 a P a i + So(»)(j)) 



5« = * 4(5 " + hi? , »;;.„ = rr;-.*,„ zr 



2J -2,? 

-"i? 1 '"ij ) /T can — ^can 



ij'TT 
can ' 



(50) 



where /i?- T and 7r*{J T are transverse traceless, e.g., /iJ T — 
fiEj = 0, and 7r*{ n is related to a vector potential V^ an by 
= YLnj + - ^V^. The gravitational con- 
straints can now be solved for "J/ and 7Tcan> leaving hj? and 
7r cari T as tne nna l degrees of freedom of the gravitational 
field. Notice that our gauge condition 7r" an = deviates 
from the original ADM one 7r" = by spin corrections at 
5PN. The action reads, with 167r-abuse of canonicity, 



W = I d 4 x^^h 



can t TT 
l07T j ' 



dt 



(00')' 



WO) _ e 



(51) 

and is in fully reduced canonical form. The dynamics is 
completely described by the ADM energy E, which turns 
into the volume integral 

E = ~J d 3 xA^[z\ P h S m) , hJ?,*W?) , (52) 

and is the total Hamiltonian (E = H), once it is expressed 
in terms of the canonical variables by solving (|4"5)l with 

The equal-time Poisson bracket relations take the stan- 
dard form (x = (x 1 )), 



{z l ,P,j} — 5ij , S^)} — £i]kS( k ) 



(53) 



{^.T (X)t)) ^TT (x ^ )} = 167r5 X™ J(x _ x0) (54) 

zero otherwise, where = §e«jk<%)(fe)j e vk = (i — j)(J~ 
k)(k — i)/2, and is the TT-projection operator, see, 

e.g., [TJ]. The Hamiltonian H generates the time evo- 
lution in the reduced matter-field phase space. Gener- 
alization and application to many-body systems is quite 
straightforward, see [12]. The total linear (P* ot ) and an- 
gular (Jfj 1 ) momenta take the forms (particle labels are 
denoted a), 



ptot _ p . \_ [ c ]3 T7r klTT h TT 

r i — / -^g^ / CI ^"can n kl,i > 



(55) 



8tt 



H 3 ( itTT, TT _ jfcTT, TTa 
a x (."can "fcj "can n ki ) 



(56) 



1 J .3 / ^ fciTT? TT j fc/TTi TT\ 
d X(X 7T can ft fcZj — X 7T can fl/j^J . 



and are obtained from the reduced action in the standard 
Noether manner. 

Correctness of the developed formalism has been ex- 
plicitly checked through conservative 3PN and dissipative 
3.5PN orders (or conservative 3.5PN and dissipative 4PN 
orders for spin-orbit interaction if spin is counted of or- 
der 1/c) by an independent method based on the full 
Einstein field equations with the Tulczyjew-Dixon stress- 
energy tensor as source term, [12(113] . Further, an alterna- 
tive derivation of P* ot and J*?* via surface integrals and 
the momentum constraint is given in [13] . which provides 
a check of the canonicity of the variables up to all orders. 

Note added in v3: Extensions of the action in [22] 
to gravitational interactions have already been considered 
in [28] . Whereas we have only given a minimal coupling 
here, Eq. ([2]), in [28] even nonminimal couplings to gravity 
and couplings to the electromagnetic field were discussed. 
However, a separation into independent variations for field 
c5e QAl and matter SA as necessary in the present paper 
was not considered. 



We thank S. Hergt for helpful discussions. This work 
is supported by the Deutsche Forschungsgemeinschaft 
(DFG) through SFB/TR7 "Gravitational Wave Astron- 
omy" and GRK 1523. 

REFERENCES 

[1] P. A. M. Dirac, in Recent Developments in General Rela- 
tivity (Pergamon Press, Oxford 1962), p. 191. 

[2] T. W. B. Kibble, J. Math. Phys. 4, 1433 (1963). 

[3] S. Deser and C. J. Isham, Phys. Rev. D 14, 2505 (1976). 

[4] J. Geheniau and M. Henneaux, Gen. Relativ. Gravit. 8, 
611 (1977). 

[5] J. E. Nelson and C. Teitelboim, Ann. Phys. (N.Y.) 116, 
86 (1978). 

[6] M. Henneaux, Phys. Rev. D 27, 986 (1983). 

[7] H. P. Kiinzle, J. Math. Phys. 13, 739 (1972). 

[8] E. Barausse, E. Racine, and A. Buonanno, Phys. Rev. D 
80, 104025 (2009). 

[9] K. Yee and M. Bander, Phys. Rev. D 48, 2797 (1993). 
[10] R. A. Porto, Proceedings of the 11th Marcel Grossmann 
Meeting on General Relativity. Part C, eds. H. Kleinert 
and R. T. Jantzen (World Scientific, New Jersey 2008), 
p. 2493; R. A. Porto and I. Z. Rothstein, Phys. Rev. D 
78, 044012 (2008); R. A. Porto and I. Z. Rothstein, Phys. 
Rev. D 78, 044013 (2008). 
[11] R. Arnowitt, S. Deser, and C. W. Misner, in Grav- 
itation: An Introduction to Current Research, ed. L. 



p-5 



J. Steinhoff and G. Schafer 



Witten (Wile y, New York 1962), p. 227 (reprinted in 
arXiv:gr-qc/0405109j ). 
[12] J. Steinhoff, G. Schafer, and S. Hergt, Phys. Rev. D 77, 
104018 (2008). 

[13] J. Steinhoff and H. Wang, Phys. Rev. D 81, 024022 
(2010). 

[14] J. Steinhoff, S. Hergt, and G. Schafer, Phys. Rev. D 78, 
101503(R) (2008); J. Steinhoff and G. Schafer, Phys. Rev. 
D 80, 088501 (2009). 

[15] B. M. Barker and R. F. O'Connell, Gen. Relativ. Gravit. 
11, 149 (1979); Phys. Rev. D 12, 329 (1975). 

[16] P. D. D'Eath, Phys. Rev. D 12, 2183 (1975); K. S. Thorne 
and J. B. Hartle, Phys. Rev. D 31, 1815 (1985). 

[17] T. Damour, P. Jaranowski, and G. Schafer, Phys. Rev. D 
77, 064032 (2008). 

[18] G. Faye, L. Blanchet, and A. Buonanno, Phys. Rev. D 74, 
104033 (2006); L. Blanchet, A. Buonanno, and G. Faye 
Phys. Rev. D 74, 104034 (2006); 75, 049903(E) (2007). 

[19] J. Isenberg and J. Nester, in General Relativity and Grav- 
itation. One Hundred Years After the Birth of Albert Ein- 
stein, Vol. 1, ed. A. Held (Plenum Press, New York 1980), 
p. 23. 

[20] T. Ledvinka, G. Schafer, and J. Bicak, Phys. Rev. Lett. 

100, 251101 (2008). 
[21] Y. N. Obukhov, A. J. Silenko, and O. V. Teryaev, Phys. 

Rev. D 80, 064044 (2009). 
[22] A. Hanson and T. Regge, Ann. Phys. (NY.) 87, 498 

(1974). 

[23] R. A. Porto, Phys. Rev. D 73, 104031 (2006). 
[24] J. Schwinger, Phys. Rev. 130, 1253 (1963). 
[25] G. N. Fleming, Phys. Rev. 137, B 188 (1965). 
[26] B. S. DeWitt, Phys. Rev. 160, 1113 (1967). 
[27] T. Regge and C. Teitelboim, Ann. Phys. (NY.) 88, 286 
(1974). 

[28] I. Bailey, W. Israel, Commun. Math. Phys. 42, 65 (1975). 



p-6 



